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ABSTRACT 


In this work we discuss and compare the methods of Poincare-von Zeipel 
and Brown-Shook for the elimination of the short period terms from aHamil- 
tonian by means of canonical transformations. We develop the general theory, 
and also supply a table of operators which serve both to eliminate short 
period effects of any order and to invert the canonical transformation, i.e., to 
express the original elements in terms of the elements affected only by the 
long period perturbations. The generating functions which produce these op- 
erators are the Taylor and the Lagrange operators. The Lagrange operator 
is a generalization of the classical operator to the case of several independent 
variables. In every attempt to solve the problem of the general perturbations 
using the electronic machines, the table of differential operators and the 
recurrence relations between them constitute an essential part of the 
programming. 

The elimination of the short period terms seems to be somewhat simpler 
in the Brown-Shook method than in Poincare-von Zeipel method. However, 
the process of determination of the original, osculating, elements in terms of 
the elements affected by the long period perturbations only is more compli- 
cated in the Brown-Shook method than in von ZeipePs method, because the 
inversion of the canonical transformation in the Brown-Shook method requires 
the application of a chain of Lagrange operators. In the Poincare-von Zeipel 
method, the use of only one Lagrange operator is required for the process of 
the inversion. In this respect the Brown-Shook method resembles the clas- 
sical method of Delaunay. Of course, if we are interested in the first order 
effects only, then both methods coincide, and the statement of Jeffreys about 
the identity of both methods remains valid. 
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ON THE POINCARE'-VON ZEIPEL AND BROWN- SHOOK 
METHODS OF THE ELIMINATION OF THE SHORT 
PERIOD TERMS FROM A HAMILTONIAN 

by 

P. Musen 

Goddard Space Flight Center 


INTRODUCTION 

In this work we discuss and compare the methods of Poincare (1892)-von Zeipel (1916) and 
Brown-Shook (1933) for eliminating the short period terms from a Hamiltonian by means of canon- 
ical transformations. We also develop a system of differential operators to be used to perform 
this elimination up to any desired order of a small parameter. After the short period terms are 
eliminated, only secular, long period and critical terms, as caused by the near commensurability 
of the mean motions of celestial bodies, are left in the transformed Hamiltonian. We have in view 
applications to planetary or lunar problems, and the mathematical as well as physical importance 
of eliminating the short period terms in these two types of problem is quite evident. The behavior 
of a planetary orbit over a very long interval of time depends mainly upon the transformed 
Hamiltonian. In connection with this topic, the recent work by Meffroy (1967) deserves 
special attention. 

In the lunar problem, the transformed Hamiltonian serves to determine the mean motion of the 
perigee and of the node. In addition, the Poincare-von Zeipel and Brown-Shook methods because 
of their concise form, represent excellent tools in understanding the effect of the propagation of the 
influence of lower-order effects into the effects of higher orders. 

The elimination of all significant short period terms, if performed by hand, can be extremely 
time consuming. It took Delaunay (1860) about twenty years to complete his lunar theory. In the 
recent work by Hori (1963), based on the application of von ZeipePs method, a considerable portion 
of Delaunay’s theory is obtained in a shorter and faster manner. The application of von ZeipePs 
theory to the motion of the artificial satellite by Brouwer (1959) belongs to the most elegant results 
obtained in recent times. 

At the present time, the idea of applying electronic machines to obtain the analytical expan- 
sion of the general perturbations is becoming very popular. We are still far from solving the 
problem of the literal expansion in a complete form by machines, but successful experiments are 
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being performed at several institutes and the future looks promising. Thus, it is quite possible 
that in the very near future we shall witness significant progress in celestial mechanics, which will 
be comparable only to the progress in numerical methods since electronic machines were first 
applied on a large scale nearly twenty years ago. 

In every attempt to solve the problem of the general perturbations using the machines, the 
table of differential operators and the recurrence relations between them will constitute the es- 
sential part of the programming. The ’'generating functions” which produce the necessary operators 
are the Taylor and Lagrange operators. The Lagrange operator is associated with a generalization 
of the classical Lagrange expansion to the case of several variables. 

In this work, we give a general theory and a table and formulas for operators which serve to 
eliminate short period effects up to any order, and to invert the canonical transformation, i.e., to 
express the original elements in terms of the elements affected by the long period perturbations 
only. 

In any particular problem, the selection of significant operators and the programming can be 
facilitated on the basis of the formalized procedure for the von Zeipel and Brown-Shook methods. 

A part of the exposition given here is from lectures on von Zeipel’s method given by the author 
at the University of Maryland. An approach different from that presented in this work was given 
by Giacaglia (1964). The recent interesting results obtained by Hori (1966) using Lie series de- 
serve to be mentioned. 


TAYLOR AND LAGRANGE OPERATORS 

The basic operators which we use in this work are 


h(.v) • V x 



a~ 1 


and 


where 


n 



a = 1 


h - (lv h 2 , •••. h n ) 


is a vector function of the position vector 

x = (xj, x 2 , •••, X n ) 


( 1 ) 


( 2 ) 


2 



and 


V = ( J_ , _jL J_\ 

X \d x l ^ x 2 ^ X n/ 

is the gradient operator relative to x. We have by definition 

n 

(h • V x ) 0(x) = h • (V x <£) = 7" 1 h a 

a= 1 

and 

n 

(V x • h )<£(x) = V x • (h 4 >) = ^ 

a= 1 

In particular, <f> can be a vector or a polyadic. 

With a given set of vectors 

gj O) = (gjl- Bj2- Bjn) ’ j = 1, 2, 3, 

we can associate two sets of operators: 

g j ’ V x “ g j“ c?x a 

a= 1 

and 

V x ’ Sj 2^ gja • 

n = 1 


(3) 

(4) 

(5) 


(6) 


(7) 


Furthermore, for the purpose of elimination of the short period terms from a Hamiltonian, we in- 
troduce the products of operators (6) and (7). By definition: 


(si • V x)® (g 2 • V x)® "■ ® (s m • V x) 


n® 

a = 1 



d m 

g m,a dx„ * * * (9x^ 

m a j a 2 


( 8 ) 


3 



m 

( v *-g 1 )©( v x-g 2 )©---©( v x-gj = n @(v *' sj 


V~ ' <9™ 

/ <9x <9x ••• dx Si, a , g 2,a '** S m ,a » (9) 

L I a i a 2 a m 12 m V f 


a . ,a„ 


a a , # * * r a - 1 , 2 , •*' n 

12 T m ' T 


Evidently these products resemble the direct Kronesker products from the theory of represen- 
tations. It is important to emphasize that the products (8) and (9) do not mean that the operators 
(6) or (7) are being applied in succession. The operators are first being multiplied together, then 
applied to the function standing to the right. For example: 


(*i -V x)® (*2 • V x)<f> ^ Ui • V x){(s 2 • V x^)} - 


and 


( V x • *,) © ( V x • S 2 )^ * ( V x-g,){ V x(g 2 ^)} 

We can write (8) and (9) also in the form 


m 

n® 


g<* * V X 


(10) 


m 

J|© ( V X • Sa) 


= v« * 


m 

n 


(id 


where the powers of V and the product g 2 > ’ * ’ g n are P°lyadics and the asterisks mean a 
complete contraction, until a scalar is obtained. In particular 



* V 

* X 

= gj 

• V x 

V x 

* gj 

it 

< 

■ «j 


The product of the identical operators appears in elementary calculus and can be written as a 
symbolic power, 


( h ‘ v x )® (h-V x )0 -••© (h-V x ) = (h • V x )" = h" * v x ™ , 
m times 


4 


i nail 1 1 ■!■■■■■ 


hi 



and 


( V x • h )® ( V x • h )@ ••• ® ( v x • h ) = ( V x • h ) m = V x m * h m ■ 
m times 


Decomposing h along the fixed axes, 

h = (h,. h 2 , h n ) ; 


we have 


( h>V x)' 



a 

(?Xj <9x 


d m 
a 2 . 


<9x 


( V x ■ h) 



a. . a 0 

^x x dx 2 


dx 


h l 1 




a + a + • * • + a - m . 

1 z m 


If h is expanded into a series in powers of a small parameter 
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then the Taylor operator 


T = exp ( h * V x ) 


E 

r, - n 


( h -Vf 

a ! 



CL- 0 


and the Lagrange operator 


A - exp ( V x • h ) 


E 

rt — C\ 


( V x-h) a 

a ! 


E 


v “ * h a 
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can be expanded into series of the form 


t s L T - • a = E A ‘ ■ (i 2 ) 

a- 0 a= 0 

where T j and A. are polynomials in h k • V x , andV x • h k (k = 1, 2, • • •, j) respectively and of order 
j relative to the small parameter. We can represent the operators T. as the sum of polynomials 
homogeneous in h s * V x ( s = 1, 2, • • •)? 


T. 

j 



where T. k is of degree k relative to the operators h s * V % . The T. -operators can be associated with 
the name of Faa de Bruno (1855). Making use of the table of Faa de Bruno operators (Duboshin, 
1950) we obtain 


T o = 1 • 

T l,l = h ! ' A ' 

T 2,l = h 2 ' A • 

T 2>2 = 5 (V A) 2 = 1 **1 2 * V X 2 . 

T 3>1 = h 3 • V X ■ 

T 3,2 = (hr V x)©( h 2' V x) = h ! h 2 * V X 2 ■ 

T 3> 3 = i ( h l • V x) 3 = F h l 3 * V X . 

T 4.1 = h 4 • A • 

T 4 , 2 = (h, • V x)® ( h 3 • V x) + i ( h 2 • V x) 2 = [\K + l h 2 2 )*V 

T 4 ,3 = T(V V x) 2 ©(V V x) = J h ! 2 h 2 * A 3 ■ 

T 4, 4 = J4 ( h l ' V x) 4 = 14 h l * V ’ 
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I 


T 


5,2 


T 


5,3 


T 


5,4 


T 


5,5 


T 


6, 


1 


T 


6,2 


T 


6,3 


T 


6,4 


T 


6,5 


T 


6 , 6 


T 


7,1 


T 


7,2 


= ( h . • V x)0 (V V x) + K- V x)©(V V x) 

= (h, h 4 + h 2 h 3 ) * V x 2 , 

= l( h 1 - V x) 2 ©( h 3- V x) + l( h !- V x) ®( h 2 ' V x) 2 

= i (h, 2 h 3 + h a h 2 ^) * V x 3 , 

= F ( h l • V x) 3 © 0*2 ' V x) = F b l h 2 . 

= 120 ( h i ' V x) 5 = 120 h i 5 * V ’ 

= h a ' V x • 

= 0*1 • V x)© ( h 5 ' V x) + ( h 2- V x)©( h 4‘ V x) + l( h 3‘ V x) 2 

= ( h l h 5 +h 2 h 4 + \ h 3 2 ) * V X 2 ' 

= l(hi-V x ) 2 ©( h 4 'V x ) + (h 1 .V x )(h 2 -V x )(h 3 .V x ) 4(h 2 -V x )3 

(2 hj 2 h 4 + hj h 2 h 3 + -g h 2 3 ) * V x 3 , 

= F ( h i • \) 3 © (1*3 • V x) + i (•*! • v x) 2 © (h 2 • v x ) 2 

= (f h l 3 h 3 + i h l 2h 2 2 ) * V X 4 > 

= ^( h l ' V x) 4 ® ( h 2 ' V x) = ^" h l 4h 2 * V X 5 ' 

= 720 ( h i ' V x) 6 = 720 h i 6 * V x 6 ' 


0*1 • V x)©(V V x) + 0> 2 - V x)©(V V x) + (V V x)® 0*4* V x) 

( h l h 6 +l *2 h 5 + h 3 h 4) * V X 2 - 
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<J\\ h- 


T 7.3 = T( h l- V x) 2 ©( h S- V x) + ( h l' V x)® (V V x)®(V V x) 

+ T ( h l ‘ V x) © ( h 3 • V x) 2 + T ( h 2 • V x) 2 © ( h 3 • V x) 
= (t h l 2 h s + H 1 h 2 h 4 + T h l h 3 2 + T h 2 2 h 3 ) * V X 3 . 

T 7.4 = F ( h l ' V x) 3 © (K ■ V x) + i (K ■ V x) 2 © ( h 2 ■ V x) © ( h 3 • V x) + i (»>, * V X ) © ( h 2 ' V X ) 3 

= (f h l 3 K + i h l 2 h 2 h 3 + F h : h 2 3 ) * V - 

T 7 .s = ^( h i- V x) 4 ©( h 3- V x) + i (V V x) 3 ©(V V x) 2 

= ("24 h 4 h 3 + 12 h i 3 h 2 2 ) * V x 5 ’ 

T 7,6 = m ( h .- V x) 5 ©( h 2 - V x) = no h l 5 h 2 * V X 6 ' 

T 7 ,7 = 5040 ( h . - V x) 7 = 5C*0 S 7 * V x . 

T 8, 1 ~ ^8 ' ’ 

T 8.2 = ( h , • V x)® 0*7 ’ V x) + { h 2 - V x)©( h 6 ‘ V x) + ( h 3 - V x)©( h 5 - V x) + l( h 4’ V x) 2 

= (S'>7 +h 2V h 3 h 5 + T h 4 2 )* V X • 

T 8,3 = l( h 1 ' V x) 2 ©( h 6- V x) + ( h l- V x)® C> 2 - V x)©( h 5- V x) + ( h l' V x)© ( h 3- V x)©( h 4- V x) 

+ F( h 2- V x) 2 ©( h 4' V x) + ?C' 2 - V x)©( h 3- V x) 2 

= (y S 2 h 6 + S h 2 h 5 + h l h 3 *>4 + \ h 2 2 h 4 + F h 2 «* 3 2 ) * V X ■ 

( h i • V x) 3 © («* 5 • V x) + \ Oh ' V x ) 2 © (h 2 • V x )©(h 4 • V x ) + i (h, • V x )2© (h, • V x )2 

+ i ( h i - V x) © Oh - V x) 2 ©( h 3 - V x) + A(»* a -V x y 


8 



(f h . 3 h s + 1 h ! 2 h 2 h 4 + ■? V h 3 2 + i h I h 2 2 h 3 + M h 2 4 ) * 


V x 4 


T 8.5 = -54 (V V x) 4 ©(V V x) + i ( h ! • V x) 3 © ( h 2 • V x)@ (V V x) + T5 (K - V x) 2 © (V V x) 3 

= (^4 h 4 h 4 + F h l 3 h 2 h 3 + TS h l 2 h 2 3 ) * V X 5 - 

T 8, 6 = 150 ( h r V x) 5 © ( h 3 ’ V x) ( h l' V x) 4 ©( h 2 - V x) 2 

= (l30 h l S|, 3 + iV h 2 2 ) * V X ' 

T 8,7 = 750 (V V x) 6 ®(V V x) = 750 h i 6 h 2 * V x 7 . 


8 , 8 


40320 ( h i ' V x) 8 


- h 8 * V 8 

40320 h i * V x 


The recent work of Sconzo and Valenzuela (1967) on the computation of Faa de Bruno operators 
on the electronic machines must be mentioned. 

In the process of elimination of the short period terms from a Hamiltonian, we shall make use 
of the following form of the Taylor expansion: if 


F(x) = 

a ~ 0 


and 


h = ha (x) 
a = 1 


then making use of (12) we obtain: 

oo a 

F(x+h) = TF(x) = 2^ 2Z T “-' 3 F “ (X) • 

a- 0 /3 - 0 

The following two important relations exist between the operators T. k and h s • V x : 


V ( h ! • V x) a ' © ( h 2 ’ V xT 2 © ' ' - © ( h p ’ V x) ap 

Z_ a.! a'. a! 


( 13 ) 


9 



where 


p 



s = 1 



j . 


kT 


j, k 


j-k+1 

E 


cr~ 1 


h 


cr 


© ^ j-cr.k-l 


(14) 


The expansion of the Lagrange operator 


A - exp (V x • h ) 


is performed in a manner similar to the expansion of T. We obtain: 


a • !>• 

a= 0 


and 


A. 

j 



k = 1 


where A. k are the homogeneous polynomials of degree k relative to the operators V x • h s . We have, 
for example: 


A 2.2 = 1 (V K) 2 - I V X 2 * W 

A 3,l = V X ' h 3 ’ 

A 3.2 = ( V x' h l)©( V X-h 2 ) = V X 2 * h ! h 2 

a 3. 3 = 5(V»«i) 3 = i v * 3 * V • 


10 



V X • h 4 - 


4,2 


( V x • h i) © (h 3 • V x) + i ( V x • h 2 ) 2 
V x 2 * ( h ! h 3 + T h 2 2 ) • 


The relations (13) and (14) exist also for A j k operators. Let us assume that x and y, defined on 
the n-dimensional manifold M, are related by the equation 


x = y - h(y) , ( 15 ) 

where h(y) is analytic on M. We assume that the transformation (15) is continuous on M and also 
non-singular, so that the representation of y as an analytic function of x is possible. Let F(y) be a 
given analytic function of y, defined in M. Under this assumption a generalized Lagrange expansion 
is valid (Stieltjes, 1885), (Good, 1960), (Dederick and Chu, 1964): 

F(y) = A{h(x)}{j[h(x)] F(x)} ' (16) 


where 


A{h(x)} exp K • h (*)} 


is the Lagrange operator and 


j[h(x)] - det [i - V x h(x)] . 

Taking (16) in the expanded form we have (Dederick and Chu, 1964) 


F(y) - 2^ v x“ * {[ h ( x )] a j[ h O)] F (x)} • 

a “ 0 


In particular we have 


h (y) - r.jr v x a * {[ h ( x )] a+1 J[ h (*)]) 

a=0 


11 



and the inversion of (15) 


y = * + v * * {[ h (*)] a+1 J [ h (*)]} 

a= 0 

= x + A[h(x)] j[h(x)] h(x) . 

An essential part in the application of Lagrange formula is the expansion of Jacobians of the 
form 


J = det {i - V x h(x)} 

in powers of a small parameter. 

The expansion of determinants of the type 


det (I - A) 

for the n-dimensional matrix A with the sufficiently small norm can be accomplished by the method 
of Leverrier (1840), based on the representation of the determinant of I - A in the form: 


det(I-A) - exp tr log (I -A) 



where 


x k 


= tr ( A k ) 


Evidently the expansion of det ( I - A) will be of the form 


det (I- A) = Y] (-l) k G k • ( 17 ) 

k= 0 


where 


( _ l) k G = T. 




*2 

2 


S) ■ L 


A.+A-+ •••+A. A, A 0 

( 1 ) Xi x 2 


(1 1 2 2 


5 )(V 


v- 


■K') 


( 18 ) 
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and where the polynomials (-l) k G k are obtained from T k simply by replacing each h. • V by -(y./j). 
The polynomials G k in Higher Algebra are associated with the name of Warring. 

In the explicit form we have 

det ( I - A) = 1 “ Yi + \ [xi ~X 2 ) “ F [x? " 3 *i * 2 + 2 V 3 ) + lA (xi ~ 6X1 X 2 + 8 Vi X3 + 3 v 2 2 “ 6 X A ) 

■ 150 (*i " 10 ^1 3 *2 + 20 ^1 2 *3 + 15 *1 *2 ~ 30 *i *4 " 20* 2 y 3 + 24y s ) + • • * . 
If A is expanded into a power series relative to a small parameter 


00 



j= 1 


then the procedure should be slightly changed. We have 


log 




Y', IT = - X 1 - ( X 2 + i X L 2 ) 

k= 1 


-( 


X 3 +X 1 


X 2 + 


I X l 3 ) ' ( X 4 + 


+ X, 


■.«) 


( 19 ) 


Taking into account that tr (MN) = tr (NM), we conclude that in order to obtain the expansion of 
tr log (I - A) we simply replace in the* last equation each product of the form 


by 




Thus we obtain 


det (I - A) - exp +X P 2 + '/' 3 + ’ * *) 




k=0 


( 20 ) 
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where we set 


0 1 x - tr Aj , 

^2 = “ tr ( A 2 + 5 A i 2 ) ■ 

^3 = " tr ( A 3 + A 1 A 2 + 1 A l 3 ) ' 

>A 4 ‘ - tr ( A 4 + A . A 3 + i A 2 2 + V A 2 + i A l 4 ) • (21) 


If necessary, the extension of expansion (19) and the table of traces given above can be easily ac- 
complished by means of the recurrence formula 


c 


n + 1 


n 

£ 


(n + 1) x n 


In the theory of the elimination of short period terms from a Hamiltonian, all matrices A j are 
trigonometric series and the diagonalization of I - A cannot be done easily. For this reason the 
Leverrier method, because it requires only the multiplication of matrices and the formation of 
traces, appears to be a very convenient way to obtain expansion (20). 

If h(x) is expanded in powers of a small parameter, 


03 

h ( x ) h a ( X) ’ 


then setting 


\ VK oo ■ 


in (21) we deduce the expansion of j{h(x)} 


j 22 j - • 

a - 0 


where J a (a = 0, 1, 2, • • •) are obtained from the formulas for the T a -operators by replacing 


14 



by 0i 9 0 2 9 * * *> respectively. If F(y) is also given as an expansion in powers of the small 
parameter, 


F (y) - 2 ^ F a (y) 

a = 0 


then the expansion of F(y) and in powers of the same small parameter has the form 


In particular 


F(y) = A{h(x)} {j [h(x)] F(x)} 


00 a 

a= 0 /3 = 0 y~ 0 


00 a yS 


y - x + 


YLYLH A “^ J /3-y h r ^ 


a=0 /3= 0 y= 0 


( 22 ) 


(23) 


ON THE METHOD OF VON ZEIPEL 

In developing von Zeipel T s method of elimination of the short period terms, we follow the work 
briefly described previously by the author (Musen, 1965). 

Let us consider the system of canonical equations 


dx 

dt = 

+ V y F , 

II 

'S'te 

- V * F 

df = 

dt 

+ F , 

d v - 

dt " 

- V £ F 


(24) 


where the vector x has the same number of components as y, and f has the same number of com- 
ponents as . 

We assume the Hamiltonian F to have the form 

F(x, f; y, 7?) = F 0 (x) + FjCx, y, tj) + F 2 (x, y , v ) + F 3 (x, y , v ) + * * ■ ’ (25) 


15 



where 


F j (*, ('< y- v) - Ara'-.y (*> i) exp i(m • y + p • 77) , ( 2 6) 

m, m 

and m and n are the vectors whose components are integers. 

The terms in (25) or (26) are classified: 

as short periodic if m / 0, 
as long periodic if m = 0, but i± / 0 
as purely secular if m = 0 and ^ = 0. 

The basic idea of the Poincare (1892) and von Zeipel (1916) methods of elimination of the short 
period terms from (25) consists of finding such a canonical transformation 


(x, y; v) - ( x *, y*; v*) . 

x - x* + V y S(x*, f*; y,v) > (27) 

y* = y + v x * s ( x *’ £*; y. v) . (28) 

f r g* + Vtj S(x* , f*; y, v) > (29) 

= 7 ? + S( x* , £*; y, 77 ) , ( 30 ) 


where S is of the order of perturbations, such that no new short period argument y* is present in the 
transformed Hamiltonian F*. Thus 


F(x, g; y, 7?) = F* (x*, f 77*) . 

Substituting (27) through (30) into the last equation we obtain 

f{x* +v y S(x* , f*; y, 77) , g* +^ v S(x*, f * ; y, 7;) ; y, 77} = F* {x*, f*; 77 + S(x*, g * ; y, 17)} ’ ( 31 ) 

which serves for the determination of S and the form of F*. In the process of solving ( 31 ), the 
notation does not play any essential role and thus the asterisks can be omitted. We write ( 31 ) in 
the form 

F {x +V y S(x, g; y, 77) , g + S(X, g- y, 77) ; y, t; } = F* {x , g ; - , 77 + S(x, g\ y , 77 )} • ( 32 ) 

In this respect we follow in von ZeipePs footsteps. 
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At this point, for the sake of brevity, it is convenient to combine the subspaces of x and f, and 
of y and 77 , into two full spaces and to introduce the combined vectors 


«) • * ■ «) 


The corresponding gradient operators are 


/ v *\ „ A 

\% ) ’ v w - yi 


Then the last equation can be re-written in the form 


where we set 


TF(u , w) = T* F* (x, 77) 


T = exp(V w S-V u ) = exp(V y S • V x + V, S • V f ) 


and 


T* - exp ( S • ) . 


Assuming the expansions 


S - Sj + S 2 + S 3 + 


F* “ F 0 * + F* + F 2 * + 


in powers of a small parameter we have 

00 00 

T = exp ^ V * Sj • V u - T j 


T* - exp Vg S. ■ V v - 2^ T j* 

J=1 j = 0 


It follows from (33) and (37) through (39): 

00 j 00 j 

■ LL t '- 


F * 
k r k 


j = 0 k = 0 


j = 0 k= 0 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 
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and thus 


Z>« f > = £ t . 

k = 0 k=0 


* F * 
-k r k 


Making use of the formulas for the expansion of the T-operators we have 


’.S, • v u 


T 2 ^S 2 • v„ + 2 ( V »S 1 -v u ) 

= V W S 2 * (V wSl ) 2 * V „ 2 


( 41 ) 


V w s 3 • y u + ( V w S i • V u )( V w S 2 • V u ) + 5 ( V w S i • V u) 3 

V wS 3 * + V W S, V W S 2 * V H > + i (W * V , 

Vw s 4 • V u + (v w s, • V u )(v w s 3 • v u ) + \ (v w s 2 • v u ) 2 + 5 (v w s, • v u ) 2 (v w s 2 • V u ) + ~ (v w s t - v u ) 4 

V w s 4 * v u + v w s 1 v w s 3 * v u 2 + \ (V w s 2 y * v u 2+ | (V wSl ) 2 V w S 2 * V u 3 + (v wSl ) 4 * v u 4 , 


In the expansions of the T. -operators, the expressions V w s k (k = 1, 2, • • •) are not affected by the 
operator v In a similar manner we obtain 


T 0 * = 1 , 

V = ^S, • V, , 


T 2 * = v f s 2 • V, + | (V f Sj -V v y 

-- VfS 2 * V, + \ (v f s ,) 2 * V , 2 , 

T 3 * = Vf s 3 • V„ + (V f Sj • V,)(v f s 2 • V,) + \ (v f Sl • V,)* 

= % s 3 * V, + Vg Sj Vg s 2 * V , 2 + i (Vg s I ) 3 * V , 3 , 
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t 4 * = ^ s 4 • v, + (v f Sl • v„)(v f s 3 • V,) + J (v f s 2 • V v y + \ (7 tf S, • V v y s 2 • V,) t £ (v f Sl • V v y 
= v f S 4 * V, + (v fSl V^S 3 ) * V, 2 + J (v,,S 2 ) 2 * V + | (v fSl ) 2 v f s 2 * V? + £ (v^S,) 4 * V„« . 


Here, s k (k = 1, 2, • • •) are not affected by the operator v v . The continuation of these tables 
does not present any difficulty. 

Let us assume that x lies in a given domain D and let 


IWj , IU2 > ^3 5 


be a sequence of vectors where components are integers. If all the numbers 


m j * V X F 0 O) • X e D 

can be considered small, then all the arguments of the form 


m j * y + • v k . 


for 


j = 1, 2, 3, 

and any n k , are called the critical ones. In planetary problems they are associated with the com- 
mensurability of the mean motions of celestial bodies. The near com men sur ability conditions 
produce the small divisors in the process of integration, and thus solving the problem in terms 
of trigonometric series with the arguments of the form m * y + n • 77 can become impossible. We can 
keep in the expansion only a finite number of the critical arguments and, normally, we do not keep 
more than one value of m j? say m 1# 

We introduce now the averaging operator M which performs the extraction from a given func- 
tion of the purely secular, long periodic and the critical terms. In addition to operator M, it is con- 
venient to introduce the operator P, which extracts only the short period terms from a given function. 
Thus, if a function $ has the form 


$ 


( x > f) exp i(m * y + 11 * 77 ) 

m , m 
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then 


m - 2^ >f jL ex P i P ' l + V ^ ex P 1 ( m j ' y * */) ■ (42) 

m j , M 

We make use of the operator P, which leaves in the functions of the form (40) the short period terms 
only, and we obtain 


PD 


L 


(X, f) exp i (m * y +/* • 77) 


where m satisfies the conditions 


(42 f ) 


m / 0 m / nij (j - 1, 2 t 3, ***) • 

From (41), and using the representation of T-operators in terms of the partial derivatives of 
S A , S 2> ' ’ m 9 we deduce: 


and 


V^Fo = 0, T.*F 0 = 0 (j = 1, 2, 3, •••) 


and consequently, 


n • V y Sl = F x F * , 

»• V 2 = i ( V y S i ’ V x) 2f o + ( T x F i -T/F,*) + (F 2 -F*) 

n • V y s 3 = [(^ySi 'V x )(V y S 2 ■ V x ) + ^ (V ySl • V x )3] Fo 

+ (T 2 Fj - T 2 * Fj* ) + (Tj F 2 - T,* F 2 * ) + (F 3 -F 3 *) 


" • V y S 4 = [(V, S, ' V x )(V y S 3 • Vj + \ (V y S 2 • V x y + £ (V y S, • V x ) a (V y S, • V x ) + £ (V y S, • V x ) 4 ] F 0 

+ (T3 Fj - T 3 * Fj* ) + (T 2 F 2 -T 2 *F 2 *) + (T,F 3 -T,*F,*) + (F 4 ~F*) , 


where we set 


n 


^* F o 
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The functions S,, S 2 , • • • must contain the short period terms only, but F*, F 2 * , • • • must 
contain only the purely secular, the long periodic and the critical terms. Thus we can set 

F* = MFj , 

“ ’ 7 , S ! ° PF 1 • 


F* = (V y S, • V^^Fo + TjF, + F 2 -T*F*} 

" • V y S 2 = Pfi ( V »S, -^^Fc + TjFj + F 2 } , 


F 3 * = M 


{[( V y S, ’ y x )(V y S 2 • V x ) + i (V y S, • V x )3] Fo 
+ T 2 Fj + T, F 2 + F 3 - (t* F * + T * F*)j , 

S 3 ^ p{[( V yS, -V x )(V y S 2 -V x ) + ^ (V y S 1 -V x ) 3 ]F 0 + (T 2 F 1 + T 1 F 2 + F 3 )} 


F 4 * = M 


n • V 


{[( V y S . ■ V x)( y y S 3 • V x ) + i (V y S 2 • V x )’ + j (v y S. • V x )> (V y S 2 • V x ) + £ (v y S, ■ ? x )<] F 0 
+ (T 3 Fj + T 2 F 2 + T, F 3 + F 4 ) - (T* F* + T* F * + T* F*)} . 

S 4 = P{[( v y S, • V x )(V y S 3 • V x ) + \ (V y S 2 • V x ) 2 + \ (V y Sj - V x )* (V y S 2 • V x ) + ^ (v y S, • V x )«] F 0 
+ (T 3 F 1 + t 2 f 2 + t i f 3 ) + f 4 J . 


where we took into consideration that the product of P and M is a zero -operator 


PM = MP = 0 


and that 


M 2 = M . 
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The problem of integration of the linear partial differential equations determining S 19 S 2 , 
evidently is reducible to the integration of the set of differential equations of the type 


n * V <b - A exp i(m*y+/x - 77) 


then in this particular case 


<f> ~ irn ex P i (m * y + A* * v) 


Returning to Equations (27) through (30), which can be written in the form 


w* w + V u *S(u*, w) , 
u = u* + V S( u * , w) , 


(43) 


we can determine u, w in terms of u*, w*; as the final result we have 


w - w* - a{w*, -V u+ S(u*, w* )}• J {-V u * S( 11* , w*)} V u+ S(u*, w*) 


and 


u - u* + a{w*, ~^ u *S(u*, w* )} j{-V u *S(u*, w* )} ^w*S( u * , w*) 


where we are in agreement with the previous notations we set 


A {w* ; -V u *S(u*, w*)} - exp{-V w * * V u+ S(u*, w*)} 


and 


j{-V u ,S(u*, w*)} - det {i + V w< ,V u ,S(u*, w*)} . 

The Lagrange operator A can also be expanded into a power series with respect to the small 
parameter. Using the system of formulas developed in the previous chapter we obtain: 

00 

a{w*, -V U *S} = ’ 

j = 0 
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where 


A 0 r 1 > 

Aj “ "V ' ^w* S i ’ 

A 2 “ V ■ V W ,S 2 + 2 (V u * ’ V w+ Si) 

= " V U* * V W* S 2 + \ ( V W* S l) 2 * V U* • 

A 3 - -V u , • V w *s 3 + (^ u * * ^ w *Sj ) (v u! „ ♦ V W]fi s 2 ) ^ (v u+ * V w+ Si ) 3 

“ V u* * ^w*^3 + ^u* * (^w*^i ^w*^) ” 6 ^u* * K*sJ 

A4 " _ ^||* ' ^W* ^4 + (^U* * ^W* ^1 ) (^U* ’ ^W* ^3 ) 


1 

+ 2 

(*».• 

V w ,s 2 ) 2 

- 1 K* 

■ V s x ) 2 

(V • 

V w. S 2) 

1 

+ 24 

(V' 7 «.s,) 4 

V,.. . * 

V,_ s „ 

+ V 2 * 

v s, v 

s. 





u* 

W* 4 

u* 

w* 1 w* 

3 





1 

+ 2 

V 2 * 

u* 

( ^W* ^2 ) 

2 - b 3 , 

2 u* 

* ( V w* S l) 

2 ( V w. 

S 2 ) + J4 

V 4 * 
u* * 

( V w.S x ) 4 , 


and where 


S j = S. (u*. w*) , ( j = 1,2,3, — ) . 

The expansion of j{-V u * s} can be obtained using (20) and (21). Setting in these equations 


Aj ^w* ^ u * S. , 

tr {( V w* V «‘ S ki ) j *(v w .V u .S k2 ) i2 ••• (V W .V„, S ks ) i3 } 

*1 = x, (1 > . 

*2 = x 2 (1) - 7 *i <2> • 


(i 1 • j 2 
*k, ,k 2 „ 


•*> k s 
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^3 = X3 (1> - Xl<V> ^Xl (3> • 

^ - x 4 (1) -x^ 1 ’ -ix 2 (2) + x# l> -ixi (4) . 


we obtain the following expansion in powers of a small parameter 

j {-V S ) = 1 + 'Pi + (^2 + i '/’l 2 ) + (^3 + 'Pl'p2 + 6 'Pi) 

+ ( 04 + </>l S^3 + T ^2 2 + 1 ^1 ^2 + 24 'Pi) 1 


Jo + Jl + J 2 + 


and the expansion of the original variables u and w in powers of the small parameter and in terms 
of the variables u*, w* in accordance with (23) becomes 


qp a p 

w = w * " 2Z HH A “-/3 J ^-v v u* s r< u *’ w *> • 

a =0 f 3 = 0 y~ 0 

co a /3 

u = u * + Ka ~ p J<3 ~ y Sy (u *’ w * ) • 


THE METHOD OF BROWN AND SHOOK 

In von Zeipel’s method, the transition from a given set of canonical variables (x, f ; y, ?/) and 
the Hamiltonian f(x, ( ; y, rj) to a new set (x*, (*; y*, q*) and the Hamiltonian F* (x*, f*; 77*), 

which are not affected by the short period effect, is accomplished by a single canonical transformation. 

In the method of Brown-Shook (1933) such a transition is accomplished by a chain of canonical 
transformations. The first transformation eliminates the short period effects of the first order 
from the Hamiltonian; the second transformation eliminates the short period effects of the second 
order, etc. We continue this process till all the significant short period effects are eliminated. 

In this chapter we discuss the first and the successive steps of the elimination procedure and 
the type of operators associated with it. 

Let the original set of canonical variables be 

( x o * € 0 ; yo * ^0 ) * 

and the corresponding Hamiltonian expanded in powers of a small parameter be 

" ^0 ( x o’ ^o ; y<r ^0) - F 00 (*o) + ^01 ( x 0’ y<r ^0) + ^02 ( x 0’ ^0 * yo f ^0) + 
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Evidently y 0 is the short period argument, and t? 0 is the long period argument. We determine 
the canonical transformation 


( x o« 

fo; yo 

■ Vo) 

- 

( x i 

■ *i 

: yi 

> 71) 

X 1 = 

x o - 

v *. 

Si 

( x o 

’ £<) 

; yi 

- *7 1 ) 

fl = 

fo - 


Si 

( x o 

f 0 ; 

y^ 

’ll) 

y 0 = 

yi ~ 

V *o 

s. 

( x o- 

f 0 ’ 

y^ 

’ll) 

^0 = 

*71 “ 


Si 

( x o- 

fo! 

y^ 

’ll) 


in such a way that the Hamiltonian takes the form 


( 44 ) 


f 0 (x 0 , fo- y 0 . vo) - F i (*i- fi; yi- vi) - Fio (*,) + F „ (*,. fi •>?i) (45) 

+ F i2 (*i- (i< yi- vi) + F i3 ( x i > ^1 = yi * ^l) + - 

with no first order short period terms present. 

Substituting (44) into (45) and changing the notation we obtain the functional equation 

F 0 [x, f; y-V x S, (x, £; y, ,) , V - Vg S, (x, £, y, t,)] 

= F , [x - V y S, (x, £; y, tj) , £ " Sj (x, £; y . 77) ; y , 77] , 

which also can be written as 

ex P (- V x Sj • V y - Sj • V„)F 0 (x. £; y, 17) ( 46 ) 

= exp (- V y S j • V x - V^ S, • Vg) F, (x, £; y, 17) , 

and which serves to determine the form of S, and F,. 


Introducing vectors 



we can write (46) as 


X ( x ) F = T ( F 
1 o r o 1 1 r 1 


(47) 


where 


x (*) 

1 o 


exp 


' V U S , 


• v -> = L 


T (!) 
1 0 , a 


= exp(-V B S-V n ) = 

a=0 


and T 0 (1 J, are Faa de Bruno operators, which in this case take a simple form: 


x (*) 

1 0,a 


_ (-i)“ 


Ks r v,r 


(48) 


x (*) 

1 1 , a 


(-i) a 

a! 


( V wS 1 ' V u )“ 


(49) 


It follows from (47) 


a 



X ( X ) F 
1 0, a~/3 r 0,J3 


L 


X (*) F 
L l,a-/3 r l,/3 


a = 0, 1, 2, 


(50) 


The first two Equations (50) give 


F 


i, o 


F o,o ( x ) 


(51) 


and 


* V y S i 


+ F, 


( 52 ) 


where 


n o " ■ V x F o,o ( x ) • 
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Introducing again the M and p operators and designating by Q the corresponding integration 
operator we obtain 


Si = QPF 0il and F lfl = MF 0fl . 
The remaining Equations (50) serve to determine 


Taking into account 


T ( 1) P = 0 

A 0,a r 0,0 U 1 


T F 

A 0,0 f 0,a 


■ 0, a 


T (!) F 

1,0 r 1, a 


we obtain from (50) for a > 1 


a+l a ~ 1 

Fo.a + ( “a!— ( V y S . ' ^Y.^ 0,0 + £ F o./J “ ) 

/3= 1 


(53) 


Fj j does not contain y. Thus the short period effects of the first order are eliminated, but f 12 , 
f 13 , • • • contain such effects. After F, ; (j = 1, 2, • • •) are determined, we can return to the 
variables Xj , f 1 , y 2 , >? t by means of the substitution 


( 


£ y v \ 

fi yi h) 


into Fj j . However, this is not actually necessary. The whole procedure of successive elimination 
of the short period effects of second and higher orders can be continued in terms of x, £, y, -q. 

The elimination of the second order short period effects from Fj is achieved by means of the 
canonical transformation: 


( X l- ^ 1 * y t - Tl) ( X 2 ’ ^ 2 ’ > 2 ' v 2 ) - 
X 2 = X 1 " V y 2 S 2 ( X l- £l - > 2 - Vi) . 

^2 “ ^1 _ ^ T 7 2 ^2 ( X l> ^ 1 ' ^ 2 ’ ^ 2 ) ’ 

= y 2 ■ v x, s 2 ( x i- y 2 - v 2 ) . 

v 1 - v 2 ~ t ( x 2 < € i > y 2 < v 2 ) ■ 
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INI 


II 


III 


with the properly chosen S 2 . We must have 

F 1 ( X 1 *^ 1 ; ^1 ’ ^ l ) ~ F 2 ( X 2 ’^ 2 ; ^2 ’ ^2 ) 

■^20 ( X 2 ) + ^21 ( X 2 * ^ 2 ’ ^ 2 ) 


+ ^2 2 ( X 2 ’ ^2 ’ “* ^2 ) + *3 ( X 2 ’^ 2 ; ^2 ’ ^ 2 ) + 


and consequently, 


T/ 2 ) Fj (u, w) = T 2 (2) F 2 (u, w) , 


where 


T, (2) = exp(-V„S 2 -V w ) = £ T i!a 

a= 1 


T 2 < 2 > - exp (-V w S 2 • V„) - T 2 ( .a ) 

and 


T (2) 

1 1, a 


(- i ) a/2 



( v „s 2 -v w ) 


a/2 


for a even 


for a odd 


T (2) 
1 2 , a 


/ \a/2 

Srr- ( v *s 2 -v u ) 

[2)' 


a/ 2 


for a even 


for a odd . 


The function S 2 is of the second order and, in fact, the development of the Taylor operators 
T 1 (2) and T 2 (2) is performed in squares of the small parameter, and the operators of Faa de Bruno 
with the odd second index can be set equal to zero. 

As before, we have 

a a 

£ T i^ F .., = £ = 0 , 1 ,- a; a = 0, 1,2, ••• . (54) 

fl = 0 p- 0 
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The first three equations give 


F l, 0 F 2 , 0 ~ F 00 ( X ) * 

F l, 1 ("• W > = F 2,l • 

"o ' ^2 + F 2.2 ” F l,2 • 


From the last equation we obtain 


Taking into account 


Equations (54) become 


Q pF i,: 


F 2 , 2 MF 1,2 


T (2) F =0 

A l, CL r 1,0 U ’ 


2, a 


a- 1 

p _ T (2) p + \ (T „ F fl -T^ 2 ^F 

M,a A l,a r 0,0 / j \ 1 1, a-/2 r 1 , (3 1 2 ,a-/3 r 2 , J3 ) 


f3=l 


a - 3, 4, 5, 


r ( 2 ) f 

1 1, a r 0,0 


( ’^T (W V x) F o.o 


(55) 


for a even and is equal to zero for a odd. The Equations (55) serve for the successive determina- 
tion of F 2 3 , F 2 4 , • • •. Hence, F 2 is obtained in the form 

F 2 = F 20 O) + F 21 ( X ’ “> 7 7) + F 22 O’ ^ “■ + F 23 O* ^ ) + *’* • 

with the first and second order short period terms removed. The general problem consists of 
removing the short period terms from 

F k ~ F k , 0 ( X k ) + F k, 1 ( X k » ^k 5 ^k ) + * * * + F k,k ( X k 1 ^k ' ~ ' ^k ) + F k,k+1 ( X k ’ ^k ’ ^k » ^k ) ’ 

k = 0, X, 2 * * * , 
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by means of the properly chosen canonical transformation 


( x k- fk- 7k- 

• ’ik) ^ 

( x k+i ’ ^k+i ’ 

>w 

’I k+l) 

x k+l = x k 

~ \+i 

S k+1 ( x k » £ k ’ 

y k +i 

- ’/k+l) 

^k+i _ ^k 

~ v, »k+i 

S k+i ( x k * ; 

y k +i 

■ V k+l) 


?k 

y k +i 

■ V X k S k+1 ( x k 

* £ k ’ y k + 1 

■ ’/k+i) 

V k 

^k+i 

’ ^ k S k+l ( X k ’ 

£k ; y k +i > 

Vk+i) 


which leads to the Hamiltonian of the form 

F k +i " F k+i (^Sc+i) + F k+i,i {\+v ^k+i ; 7 k+i) + * ” +F k+i,k+i ( x k +r ^k+i* ^k+i) + F k+i,k+2 (*k+i’ ^k+r y k+ r ^k+i) + ’ 

expanded in powers of the small parameter. The condition 


T k (k+1) F k O’ w > 


= T (k+1) p / u 

A k f k+lV U - 


W) 


must be satisfied, where 


T k (k+1) = exp(-V u S k+1 -V w ) = V" T™ 


T (k + 1) ~ exo (-V S + V ) - y T (k + 1) 

A k + 1 eX P V W°k+1 U/ / x k+l,g 

a~ 1 


and 


T (k+1) 
1 k, a 



(^U S k+1 


V w)’ 


for a = 0 (mod k) 


lo 


for a 4 0 (mod k) 


T 


(k + 1) 

k+1, a 



k+1 


V u) 


for a = 0 (mod k) 


lo 


for a ^ 0 (mod k) . 
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We deduce as before 


F k+l t a F k,a 


a = 0, 1, 2, • * • k 


n • V S + F = F 

o W °k+l ^k + l.k + l ^k.k + l 


and 


"0 * S k+1 PF k,k + l 


S k+1 = Q PF k 


F = M F 

r k + l,k+l rax k,k+l 


We continue this process till all the significant short period effects are eliminated. 


CONCLUSION 

The process of determining S 1? S a , • • • in the Brown-Shook method seems to be somewhat 
simpler than the corresponding process in von ZeipeTs method. 

If we are pursuing only the elimination of short period terms and the determination of long 
period effects, but not the determination of the original elements in terms of the final elements, 
then the method of Brown and Shook can provide very good service. However, the process of 
determination of the original, osculating, elements in terms of the elements affected by the long 
period perturbations only is more complicated than von ZeipeTs method. 

It is so because performing a chain of canonical transformations will require the application 
of a chain of Lagrange operators to accomplish the return to the original elements. The same 
chain of different Lagrange operators is to be used in the Brown-Shook method if one wants to 
form the canonical transformation from the original Hamiltonian to the Hamiltonian without the 
short period terms. In this respect the Brown-Shook method resembles the classical method of 
Delaunay. Of course, if we are interested in the first order effects only, then both methods coincide. 

In this case the statement of Jeffreys (1961) about the identity of both methods remains valid. 

If all the angular canonical elements, of y and q type, are eliminated, then the transformed 
Hamiltonian will contain only the elements of x and f type. These remaining elements become 
constants and the y, q elements become the linear functions of time. Such a result, for example, 
can be obtained in the case of an artificial satellite, because in the transformed Hamiltonian the 
order of the long period terms is higher than the order of the significant secular term. 
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Unfortunately, with the present day mathematical means, such a goal cannot be accomplished 
in all problems of celestial mechanics. In such a case, we suggest numerical integration with the 
large integration step of the canonical equations associated with the transformed Hamiltonian. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Green be It, Maryland, October 2, 1967 
124-12-03-01-51 


REFERENCES 

1. Brouwer, D., Solution of the problem of artificial satellite theory without drag. Astr. J., 64, 

1, p. 378, 1959. 

2. Brown E. W . and Shook C. A., Planetary theory p. 143-173, Cambridge, 1933. 

3. Dedrick K. G. and Chu E. L., On N -dimensional Taylor and Lagrange operators. Arch. Rat. 
Mech. and Analysis, 16, pp. 385-401, 1964. 

4. Delaunay , Theorie du mouvement de la Lune. Paris vol. I, 1860; vol II, 1867. 

5. Duboshin G. N., On the computation of the higher derivatives (russ.) Publ. Sternberg Astr. 
Inst., 18, No. 2, pp. 11-32, 1950. 

6. Faa de B runo, Note sur une nouvelle formule de calcul differential. Quarterly Journ. of Math., 
1, p. 359, 1855. 

7. Giacaglia, G., Notes on von Zeipel T s method, GSFC Document X-547-64-161, 1964. 

8. Good, I. J. , Generalizations to several variables of Lagrange expansion, Proc. Cambridge 
Phil. Soc., 56, pp. 367-380, 1960. 

9. Jeffreys H ., On two methods of Brown and Shook, M. N. Royal Astr. Soc., 122 , pp. 335-338, 
1961. 

10. Hori, G. , A new approach to the solution of the main problem of the lunar theory, Astr. J., 68, 
3, pp. 125-146, 1963. 

11. Hori, G ., Theory of general perturbations with unspecified canonical variables, Jap. Astr. 
Journ. 18, No. 4, pp. 287-296, 1966. 

12. Leverrier , U. J. J., Sur les variations seculaires des elements elliptiques des sept planetes 
principales, J. de Math. (1)5, pp. 220-254, 1840. 

13. Meffroy, J., On the elimination of the short period terms of a first order general planetary 
theory through von ZeipePs method. GSFC Document X-641-67-318, 1967. 


32 



14. Poincar e, H. , Les methodes nouvelles de la mecanique celeste, vol. n, p. 17, Paris, 1892. 

15. Sc onzo, P. and Valenz uela D. High order time derivatives of powers of the radius vector. 
IBM report. Cambridge. 1967. 

16. Stieltj es, T. J., Sur une generalisation de la serie de Lagrange, Ann. Sc. Ec. Norm. Paris (3) 
2, pp. 93-98, 1885. Oeuvres, 1, pp. 445-450. 

17. von Z eipel, H., Ark. Astr. Math. Fys. 11 , No. 1, 1916 


NASA-Langley, 1968 19 


33 



FIRST CLASS MAIL 


National Aeronautics and Space Administration 

WASHINGTON, D. C. 
OFFICIAL BUSINESS 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


• i *' 1 'i 

1 r ■ ‘ • C 
i i i •: i > 


ijf * 3 ^ OBOhA 00905 

_ . ** : , s L r x “» : '•*•1 A I 0 'X Y / A r W L / 

L? ,r w ^CXICU <?7L 1 < 


I * \ / ■«' l I . / 


T r CH 4 i 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


" The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 

TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 
under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS; Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C 20546 



